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ABSTRACT 

Infrared photometry of the probable triple WC4(+0?)+08I: Wolf-Rayet system 
HD 36402 (= BAT99-38) in the Large Magellanic Cloud (LMC) shows emission char- 
acteristic of heated dust. The dust emission is variable on a time-scale of years, with 
a period near 4.7 yr, possibly associated with orbital motion of the 08 supergiant and 
the inner P ~ 3.03-d WC4-I-0 binary. The phase of maximum dust emission is close 
to that of the X-ray minimum, consistent with both processes being tied to colliding 
wind effects in an elliptical binary orbit. It is evident that Wolf-Rayet dust formation 
occurs also in metal-poor environments. 
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1 INTRODUCTION 

Population I Wolf-Rayet (WR) stars lose mass at prodigious 
rates (> lO~^M0yr~^) in dense, fast stellar winds giving the 
stars their characteristic emission-line spectra: they repre- 
sent the last stable stage in the evolution of massive stars. 
Their high luminosities and conspicuous spectra facilitate 
identification in external galaxies, while the stellar winds 
have a profound influence on the local interstellar medium, 
enriching it with nuclear-processed material and frequently 
forming bubbles. When WR stars are members of binary 
systems with other luminous stars possessing stellar winds, 
the collision of the two stellar winds can result in shock- 
heating of stellar wind material and particle acceleration, 
observable in X-ray and non-thermal radio emission. Some 
colliding-wind binaries (CWBs) incorporating WC type WR 
stars also form carbon dust in their winds, observable by its 
conspicuous infrared (IR) emission. These phenomena often 
vary periodically, especially if the stars move in an elliptic 
orbit and their separation varies so that the wind-collision 
region (WCR) moves periodically between denser or rarer 
regions of the stellar winds. Also, depending on the incli- 
nation of the system, our line of sight to the WCR passes 
through varying amounts of stellar wind as the orbit pro- 
gresses, causing variations in the free- free and photoioniza- 
tion extinction to the observed radio and X-ray emission. 
The prototype is WR140 (= HD 193793), whose IR, radio 
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and X-ray variations were linked to the motion of its WC7 
and 05 components in their 7.94-yr orbit by Williams et al. 
(1990). This system has continued to attract observations in 
all wavelength domains - for a recent summary, see Williams 
(2011a) - as a laboratory of particle acceleration (e.g. Pit- 
tard & Dougherty 2006) and led to the identification of other 
variable dust-making WR CWBs in the Galaxy. The obser- 
vation of rotating dust 'pinwheels' around some of the WR 
systems which show persistent dust emission (e.g. WR 104, 
by Tuthill et al. 2008) associates these with CWBs too. The 
persistent dust makers are nearly all WC9 stars and their 
locations show a strong preference for the metal-rich region 
of our Galaxy, generally towards the Galactic Centre. On 
the other hand, the variable and episodic WR dust mak- 
ers have spectral types ranging from WC5 to WC8 and a 
less concentrated galactic distribution - but all of them are 
located within the Solar circle. This suggests a strong in- 
fluence of environmental metallicity on the phenomenon of 
dust formation by WR systems. 

In the lower metallicity environments of the Magellanic 
Clouds, the fractions of WR/0 stars and distributions of 
WR spectral subtypes are very different from those in our 
Galaxy. In particular, of the 24 WC/WO stars in the LMC, 
23 are of type WC4 and one a W03 (Moffat 2008). The 
question arises: do any of these make carbon dust? This 
can be tackled by a survey of their spectral energy distribu- 
tions (SEDs) or a search for variability in the IR character- 
istic of the variable or episodic dust-makers like WR 140. To 
begin, we examined the Ks-hanA photometry of the LMC 
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WC4 stars measured in the Deep Near Infrared Survey of 
the Southern Sky (DENIS; Epchtein et al. 1999, Cioni et al. 
2002), Two Micron All Sky Survey (2MASS, Skrutskie et al. 
2006) including the 6X w/LMC/SMC Point Source Working 
Database, and InfraRed Survey Facility (IRSF) Magellanic 
Cloud (Kato et al. 2007) surveys for evidence of variabil- 
ity. Of the 21 LMC WC4 stars with at least three obser- 
vations, the best candidate for variability, after discounting 
stars with uncertain data or lying in very crowded regions 
where close companions could be included in the apertures 
of some but not other surveys, was HD 36402 {= BAT99- 
38, Breysacher, Azzopardi & Tester 1999). It would cer- 
tainly be worth repeating the search for variability when the 
VISTA MageUanic Clouds Survey (Cioni et al. 2011), with 
its higher resolution and homogeneity, is complete. The case 
for HD 36402 being a WR dust maker was strengthened by 
its mid-IR (3.6-8 ^m) photometry in the SPITZER SAGE 
LMC survey (Bonanos et al. 2009), where it was identified as 
the reddest WR star. A preliminary account of the present 
study was presented at the 2010 LIAC (Wilhams 2011b, Pa- 
per I), where it was shown that the 2-8 \ira SED was char- 
acteristic of of WR dust, but that at wavelengths beyond 
10 \im, the emission from the WR star could not easily be 
resolved from that by two nearby luminous young stellar ob- 
jects (YSOs) discovered by Chu et al. (2005) and that the 
70- and 160-vi.m fluxes tentatively ascribed to a ring nebula 
associated with HD 36402 by Bonanos et al. came entirely 
from the YSOs. 

Independently, Boyer et al. (2010) came to a similar 
conclusion in their study of the far-IR emission from three 
dust-making evolved stars in the LMC with Herschel instru- 
mentation. They presented the SEDs of HD 36402 and of a 
strong far-IR source, HD 36402 IRl immediately to the west 
of the WR star, having Td ~ 64 and 8 K. They considered 
the possibility that IRl could be heated by the stellar flux 
from HD 36402 but concluded that it was more likely to be 
the YSOs embedded in a molecular cloud. 

The field is certainly crowded at longer wavelengths, 
making interpretation of long- wavelength observations with 
smaller instruments (e.g. IRAS or MSX) with their larger 
beam sizes difficult, but there is no confusion in the IRAC 
bands - even at 8 (im, HD 36402 is 1.7 mag. brighter than 
the closer YSO (YSO-2, 7 arcsec to the W of the WR star). 
At this wavelength, the IRAC beam (Fazio et al. 2004) is 
smaller than 2 arcsec full width half maximum (FWHM), so 
we consider YSO-2 to be well separated from HD 36402 in 
all IRAC bands as the magnitude differences are greater at 
shorter wavelengths. The field near HD 36402 at 3.6 (im is 
shown in Fig.[T] 

The spectrum and orbit of HD 36402 were studied by 
Moffat, Niemela & Marraco (1990), who derived an orbit 
from the WR emission lines having a period of 3.03269 d 
and radial velocity (RV) amplitude K = 275 ± 17 km s~^. 
The absorption-line RV amplitude was not significantly dif- 
ferent from zero {K = 8 ± 9 km s~^), suggesting that the 
08 supergaint did not participate in the 3-d orbit and that 
HD 36402 was a triple system like 6 Mus (= WR48) or 
R 130 in the LMC. From five high-resolution lUE spectra, 
Konigsberger, Moffat & Auer (2003) measured RVs of photo- 
spheric absorption lines which fitted a sinusoidal variation in 
anti-phase with the emission-line spectrum, suggesting for- 
mation in the close companion to the WC star, but more 



Table 1. Near-infrared photometry of HD 36402 



Survey or source 


Date 


J 


Ks 


aKs 


DENIS strip 4963 


1996.92 


11.63 


10.55 


0.02 


DENIS shifted 




11.73 


10.69 


0.02 


2MASS all-sky 


2000.08 


11.76 


10.97 


0.02 


2MASS 6X PSWDB 


2000.98 


11.73 


10.82 


0.02 


2MASS 6X PSWDB 


2001.09 


11.69 


10.73 


0.02 


IRSF LMC0525-6720A 


2002.88 


11.80 


11.22 


0.01 


N51D study (ISPI) 


2006.85 


11.73 


10.75 


0.05 



probably produced by the superposition of the stationary 
absorption lines from the 08 supergiant and rotationally 
broadened absorption lines produced in the close compan- 
ion. In their study of HD 193793 (WR 140), Fitzpatrick, Sav- 
age & Sitko (1982) compared its /[/iJ spectrum with that of 
HD 36402. The A1550 C iv profile in HD 36402 showed nar- 
row interstellar lines from both the Milky Way and the LMC 
and two sets of circumstellar absorption lines. The broader 
set was shifted —3100 km s~^ relative to the LMC inter- 
stellar lines and the narrower pair by —1300 km . The 
same features, at the same velocities and strengths, were 
seen in another lUE spectrum taken three months earlier, 
ruling out their formation in a transient event. Niedzielski 
& Skorzyhski (2002) used the lUE high-resolution spectra 
to determine a terminal velocity of 3072 ± 160 km s~^, com- 
parable to those of other WC4 stars in the LMC and our 
Galaxy. It seems reasonable to attribute the —1300 km s~^ 
absorption to the 08I component; it is in line with the range 
of terminal velocities (955-2186 km s~^) measured for 08 
supergaints by Prinja, Barlow & Howarth (1990). 

In this paper, we present new observations of HD 36402 
from images taken for the N51D study by Chu et al. and 
post-cryo observations extracted from the Spttzer Heritage 
Archive. Recently, additional mid-IR observations have be- 
come available from the AKARI LMC (Kato et al. 2012) 
and WISE (Wright et al. 2010) surveys. We use these data 
to re-examine the IR variability of HD 36402 and derive a 
period, and consider its status as a dust-making WR CWB. 



2 OBSERVATIONS 

2.1 IR Photometry of HD 36402 from Surveys 

Near-IR photometry and dates of observation of HD 36402 
extracted from the DENIS, 2MASS and IRSF Magellanic 
Cloud surveys are collected in Table [T] and plotted in Fig.[2l 
For an object having the J — Ks colour of HD 36402, the 
difference between the IRSF and 2MASS Ks magnitudes is 
not expected to exceed 0.02 mag. (Kato et al. 2007, table 
10). The differences between DENIS and 2MASS magnitude 
scales are, however, greater. Delmotte et al. (2001) measured 
systematic shifts between the J and Ks scales of the two 
catalogues, finding means 5 = -0.10 (J) and -0.14 {Ks) in 
the sense DENIS minus 2MASS, but with significant varia- 
tion from DENIS strip to strip. These differences are applied 
to the data in Table [T] and Fig.0 to align the DENIS and 
2MASS magnitudes. 

Mid-IR photometry of HD 36402 is available from 
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Table 2. Mid-IR photometry of HD 36402 with Spitzer IRAC from the N51D study, the SAGE IRAC Epoch 1 and Epoch 2 Catalog 
and archive; the AKARI LMC Survey Catalogue and the WISE All-Sky and post-cryo surveys. The AKARI 511 and WISE VFSmpro 
ll-|j.m magnitudes are given in the same column. 



Date 


N3 


VFlmpro 


[3.6] 


[4.5] 


iy2mpro 


[5.8] 


57 


[8.0] 


S11/W3 


Source 




(3.2 |xm) 


(3.35 |J.m) 


(3.55 |j.m) 


(4.44 |a.m) 


(4.60 \im) 


(5.58 yim) 


(7.0 |xm) 


(7.76 ^m) 


(11 urn) 




2004.96 






9.71±0.05 


8.99±0.04 




8.54±0.04 




7.91±0.03 




N51D 


2005.55 






9.30±0.05 


8.70±0.04 




8.21±0.03 




7.81±0.02 




SAGEl 


2005.82 






9.18±0.05 


8.50±0.04 




8.11±0.03 




7.67±0.03 




SAGE2 


2007.42 


9.68±0.04 












7.70±0.02 




7.11±0.02 


LMCC 


2009.77 








8.96±0.05 












IRAC 


2009.91 






9.51±0.05 














IRAC 


2010.41 




9.44±0.02 






8.57±0.02 








7.15±0.02 


WISE 


2010.91 




9.20±0.03 






8.31±0.02 










WISE 



m 




Figure 1. Field of HD 36402 (marked 'WR') from the 3.6-n.m 
IRAC image, with North at the top. East to the left. The YSOs 
discussed by Chu et al. are marked 'YSOl' and 'YS02'. The sep- 
aration of YSO-2 from the WR star is 7 arcsec and it is 4.3 mag. 
fainter than the WR star at 3.6 [xm and 1.7 mag. fainter at 8.0 |J.m. 



the Spitzer, AKARI and WISE missions. Magnitudes of 
HD 36402 in the Spitzer SAGE survey of the Magel- 
lanic Clouds (Meixner et al. 2006) extracted from the 
NASA/IPAC Infrared Science Archive are given in Table [2] 
and the [3.6] and [4.5] magnitudes are plotted against date 
of observation in Fig. (2] 

HD 36402 was observed in the AKARI LMC Survey 
(Kato et al. 2012) in the iV3 (A 3.2-|am), S7 (A 7.0-^m) 
and Sn (A 11-^m) bands as AKARI-LMCC J052603.96- 
672956.7 on 2007 June 2-3. The data are included in Table 

m 

Owing to its high ecliptic latitude (87°), HD 36402 was 
observed very frequently by WISE in both the All-sky and 
Post-Cryo surveys. It was flagged as having a high proba- 
bility of variability in the All-Sky survey, so we examined 
this using the 173 'A'-quality Wl (profile-fitted, 3.4-jxm) 
magnitudes extracted from the WISE All-Sky Single Expo- 
sure (Lib) Source Table. The data covered 21 days (MJD 
55335-356 in 2010 May- June) and we began by searching 
for periodicity using the CLEAN algorithm (Roberts, Lehar 
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Figure 2. Synoptic photometry of HD 36402 from different 
sources: J and Ks magnitudes from the DENIS (marked •), 
2MASS (o) and IRSF {<» surveys, and the ISPI frames observed 
for the N51D study (©); [3.6] and [4.5] from IRAC (e), and Wl 
and W2 from WISE (*). 

& Dreher 1987) as implemented in the StarhnlQ period 
package (Dhillon, Privett & Duffey 2001). We found a peak 
at 5.99±0.04 d but the false alarm probability (FAP) calcu- 
lated from 250 random permutations of the data was 0.056 
(cr = 0.015), so the period is only marginally significant. 
We repeated the period search on the 202 Wl magnitudes 
observed during MJD 55522-539 (2010 November - 2011 
January from the WISE Preliminary Post-Cryo Single Ex- 
posure (Lib) Source Table, but found no significant period, 
so we discount that suggested by the All-Sky data. We also 
examined the data for a trend and found the Wl magni- 
tudes to be brightening at 0.0023±0.0007 mag d~^ during 

^ The Starlink software suite is currently maintained and dis- 
tributed by the Joint Astronomy Centre, Hilo, Hawaii. 
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the 21 days of observations in the AU-Sky survey. This is 
consistent with the brightening in M^l and 14^2 magnitudes 
between the All-Sky and post-cryo surveys (Table [5] and 
Fig.[2l), where the dates used are the average dates derived 
from the Source Tables. These show Wl brightening at an 
average rate of 0.0013±0.0002 mag d~^, comparable to the 
rates of the IRAC 3.6- and 4.5-jxm fluxes. 

At ll-^m the beams of the AKARI Sll and WISE VF3 
magnitudes are 5.5 and 6.5 arcsec FWHM respectively and 
we hestitated to use these data for HD 36402 in case of 
possible contamination from YS02, whose spectrum is rising 
in these relatively wide (4.6 and 8 ii.m respectively) filter 
pass-bands while that of HD 36402 is expected to be falling. 
But, as seen from the consistency of these fluxes with the 
SEDs fitted to the shorter wavelength data in Section 13.21 
below, most if not all of the flux measured in the 5*11 and 
VK3 bands can be attributed to the WR star. 



2.2 New photometry of HD 36402 

In their observation of the N51D region, Chu et al. (2005) 
drew attention to the bright image of HD 36402 close to 
YSOs they reported. As the IRAC observations for this 
study were taken on 2004 December 12, before the SAGE 
observations, we measured the flux from HD 36402 in the 
IRAC images to extend knowledge of its variability. The re- 
sults are included in Tabled 

The region was observed in the near-lR with the Cerro 
Tololo Inter- American Observatory (CTIO) Blanco 4.0-m 
Infrared Side Port Imager (ISPl) on 2006 November 7. Three 
fields including HD 36402 were observed and the photomet- 
ric calibration was tied to 2MASS observations of stars in 
the fields. The mean J and Ks magntudes are included in 
Table [1] 

Thirdly, the Spitzer Heritage Archive was searched for 
images of fields including HD 36402. Twelve observations 
in [4.5] during 2009 October 3-15 and eleven in [3.6] be- 
tween 2009 November 16 and December 14 were retrieved 
and magnitudes of HD 36402 measured. The WR star was 
generally near the edge of the field, very close to the edge 
in the case of the [3.6] observations, making sky estimation 
difficult. Both bands show brightening during the observa- 
tion sequences but not at statistically significant rates. The 
average magnitude and date for each of the two runs is given 
in Tableland plotted in Fig.H 



3 RESULTS 

3.1 Light curves and periodicity 

It is evident from the IR observations plotted in Fig[2]that 
the amplitude of variation increases with wavelength. This 
is characteristic of varying circumstellar dust emission by 
WR stars, which peaks in the 3-4 fim region, and reflects 
the increasing contribution of the dust emission with wave- 
length, which can be seen in the SED (FigUJ Paper 1, and 
Boyer et al.). 

Variation in the dust emission can arise from variation 
in the stellar flux heating the dust or in the amount of cir- 
cumstellar dust. We observe very little variation in J (ct = 
0.04 mag., Fig[2} because the dust contributes little to the 
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Figure 3. Photometry of HD 36402 phased to a period of 4.7 yr 
with zero phase set to the epoch of the DENIS observation. For 
plotting against phase, [3.6] and [4.5] magnitudes have been in- 
terpolated from the WISE Wl and W2 and the AKARI LMCC 
N3 and 57 magnitudes as described in the text. Points from the 
AKARI data are marked A, while the sources of the other data 
are as in Fig. [2] 



flux at this wavelength. The near constancy also suggests 
that the stellar radiation dominating the J-band flux is not 
significantly variable. There is some dispersion in the visible 
photometry of HD 36402 - e.g. V = llAA (Feitzinger & Is- 
serstedt 1983), 11.618±0.098 (Zaritsky et al. 2004) - but a 
concerted search for variability by Seggewiss, Moffat & La- 
montagne (1991) showed no variation greater than the pho- 
tometric uncertainties. We therefore deduce that the varia- 
tion in the IR comes mostly or wholly from changes in the 
amount of dust in the wind rather than the stellar radiation 
heating it. The dust forms in the shock-compressed stellar 
wind and inherits its velocity, moving away from the stars as 
seen in the pinwheels made by systems like WR 104, whose 
winds are relatively slow, or the expanding splash made by 
WR140 (Williams et al. 2009), whose wind is 2-3 times 
faster. As the dust moves further from the stars, the stellar 
radiation heating it is progressively diluted and the grains 
cool, so that their emission fades. Therefore, if a system is 
observed to show constant IR emission, the dissipating dust 
must be constantly being replenished by the condensation of 
new dust grains - the observed IR SED is a consequence of 
the balance of these two processes. The growth of IR emis- 
sion by HD 36402 in 2004-05 and 2009-10 indicates that the 
rate of dust formation must have been increasing for over a 
year in each case. Between these events, the rate of conden- 
sation must have fallen below the rate needed to replenish 
the dissipation to account for the fade in emission between 
the two rises. 

The similarity of the rates of brightening in 2004-05 and 
2009-10 suggests that the variation is periodic. To facilitate 
comparison of the WISE and IRAC data, we interpolated 
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[3.6] and [4.5] magnitudes from the WISEWl (3.4-^m) and 
W2 (4.6-ti.m) magnitudes by fitting black bodies to the Wl 
and W2 fluxes, calculating fluxes at 3.6 and 4.5 \im, and 
converting them to magnitudes. Comparison of these with 
the IRAC observations shows an interval of 4.7±0.2 yr be- 
tween the two rises. The uncertainty was estimated from 
inspection of trial flts, allowing for the spreads in epochs 
of the WISE observations, and is probably conservative. In 
Fig. [31 we show the Ks, [3.6] and [4.5] magnitudes phased 
to the period of 4.7 yr, with zero phase arbitrarily set to 
the epoch of the DENIS observation. The 2000-01 2MASS 
Ks photometry fitted to this ephemeris shows a brightening 
towards a possible maximum near the time of the DENIS 
observation parallelling that seen in the mid-IR, supporting 
the interpretation of the 4.7-yr shift as a period in the vari- 
ations as the 2MASS data come from an earlier cycle than 
the mid-IR data. 

The IRAC and WISE data can also be fit by the first 
harmonic (2.35 yr) of this period, but the [3.6] and [4.5] 
magnitudes determined from interpolating the AKARI N3 
and S7 observations by fitting a black body to them do not 
fit, so we discount this possibility. Higher harmomics are 
ruled out by introduction of scatter in the phased near-IR 
data or sequence of the IRAC data, so we conclude that the 
variations are best fit with a period near 4.7 yr. 

For the present discussion, we adopt a preliminary 
ephemeris for the near-IR maximum of 

T{Y) = 1996.9 4. 7Af 

Further near-IR photometry is needed is define the light 
curves, including fading, the minimum and especially the 
rise to the next maximum expected in 2015.7 to confirm the 
period. 

The phase of the AKARI-derrved [3.6] and [4.5] magni- 
tudes, also plotted in FigO places them on the falling branch 
of the light curve. Also, the phase of the 2006 ISPI Ks 
magnitude fits the fade to the faint Ks observed by IRSF. 
Whether the IRSF point represents the true minimum is not 
certain: the IRSF J-Ks colour (0.58) is close to the mean 
J-Ks (0.55±0.07) observed for other LMC WC4 systems, 
so there may be no dust contribution at that phase. This 
would imply that condensation of new dust grains would 
have ceased some time earlier to have allowed the existing 
dust to have dissipated and its emission to have faded suffi- 
ciently. 

The light curves are not yet well defined but it is ev- 
ident that the rise to maximum takes at least 1.5 yr. In 
this respect, HD 36402 more closely resembles the Galactic 
WC7-F09 dust maker, WR 137 (HD 192641), whose IR fiux 
takes about 3 yr to rise to maximum (Williams et al. 2001), 
than the prototype, WR 140, which shows a much faster rise 
(less than two months in K, Williams et al. 2009). Better 
definition of the IR maximum from further observations will 
show if it occurs later at longer wavelengths owing to the ef- 
fects of grain growth after condensation (cf. WR140). 



3.2 Spectral energy distributions (SEDs) 

To study the SED of HD 36402, we need observations over a 
wavelength interval wide enough to cover the fiux maximum, 
essentially contemporaneous near- and mid-IR photometry. 



We do not have such data, but note that the interval be- 
tween the 2005.82 IRAC and 2001.09 2MASS observations 
was 4.73 yr, and that between the 2010.41 and the same 
2MASS observations was 9.32 yr, close to our proposed pe- 
riod or double it. Assuming that the SED repeated with 
this period, we formed a SED combining these data to rep- 
resent the system near phase 0.9 and present it in Fig.|4l 
For a 4.7-yr period, the two sets of mid-IR data straddle the 
near-IR data; when the period is revised with more near-IR 
data, it should be possible to select more closely matched 
mid- and near-IR data. Owing to the breadth of the W3 fil- 
ter, the colour correction for calibration of monochromatic 
flux from the in-band magnitude can be significant and we 
adopted that for a 800-K black body from Wright et al. 
(2010, table 1); those for the other filters are less than 2 per 
cent. 

For the shorter wavelengths, we used the narrow-band 
V observed by Smith (1968), which avoids the strong WC 
emission lines falling in the V band, and I from the DENIS 
survey. Smith also derived a colour excess of Eb-v = 0.03, 
which seems low considering the location of HD 36402, but 
is consistent with the reddening, E{B - V) = 0.00 ± 0.03, 
derived by Morris et al. (1993) from nulling the 2200- A in- 
terstellar absorption feature. We de-reddened the photom- 
etry using Smith's Eb-v, equivalent to Ay = 0.11, and a 
standard reddening law, but the effects are too small to in- 
fluence the results signiflcantly. For the stellar plus wind 
continuum, we noted that the SED of the LMC WC4-I-05-6 
system HD 35517 (BAT99-28) defined by the 2MASS, IRAC 
and WISE data has a power-law form, with no evidence of 
dust emision. We shifted this to fit the short-wavelength 
photometry to provide the wind SED shown in Fig.[31 The 
difference between the phase 0.9 and wind SEDs can be fit 
by a black body - not as a physical model, but as a sim- 
ple characterisation of the dust emission. A least-squares fit 
gave Tbb = 980 ± 20 K, with r.m.s. residuals of 0.07 mag. 
The error would be greater if the near- and mid-IR photom- 
etry did not match, i.e. if the system does not repeat with 
a period of 4.7 yr. 

As noted above, we were hesitant to include the ll-|i.m 
data in case of contamination from YS02. The modelling 
was therefore made with and without the W3 magnitude. 
The difference in fitted Tbb was less than 20 K, showing 
that the 14^3 datum is consistent with the SED determined 
at shorter wavelengths and that most, if not all, of the mea- 
sured W3 flux comes from the WR star. 

To examine the maximum dust emission observed, we 
formed a composite SED from the DENIS (shifted as above) 
and 2010.91 WISE observations, which were taken 2.98 cy- 
cles apart. This lies a little above the SED in Fig.[31 but is 
omitted for clarity. From it we measure the dust luminosity 
to be 8 X 10"^'* W m~^. We can compare this with the total 
luminosity of the three stars heating the dust. For the 081 
star we adopt log L/Lq = 5.68 following Martins, Schaerer 
& Hillier (2005). We note the unseen companion to the WC4 
star in the 3-d orbit would have a luminosity comparable to 
that of the 08 supergiant if it were an 04 main-sequence 
star, but would be ~ 0.8 mag fainter in the visible - thereby 
accounting for its not being observed in the spectra taken 
for RV study by Moffat et al. For the WC4 star, we assume 
the mean luminosity from the study of LMC WC4 stars by 
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Crowther et al. (2002): logL/L© = 5.58. This gives a system 
luminosity log L/L© = 6.13. 

Assuming 48 kpc for the distance of the LMC, the dust 
luminosity is then only 0.004 of the stellar luminosity. Con- 
sequently, the emission is optically thin and the grain tem- 
perature is determined by its thermal equilibrium in the 
stellar radiation field. For grains having optical properties 
of amorphous carbon grains, the equilibrium temperature 
falls off with distance, r, from the stars as Tg cx r""'^* (e.g. 
Williams et al. 2009). In the optically thin case, the flux F\ 
observed at our distance d from the system is simply pro- 
portional to the products of the grain opacity and Planck 
functions, Hi\B[\,T^{r)), weighted according to the radial 
mass distribution mg(r). The last is determined by the past 
history of dust formation: in the simplest case of constant 
dust formation, it would take the form mg(r) oc , but in 
the case of our SEDs on the rise to maximum, the gradient 
would be steeper, with relatively more hot dust nearer the 
stars. We therefore defer modelling the SEDs until we have 
a more complete picture of the dust formation history from 
near-IR data, together with sets of mid-IR and near-IR data 
more certainly at the same phase. 

To sample the fading branch of the light curve, we 
formed a SED from the AKARI iV3, S7 and 511 data, ob- 
served near phase 0.2 on our ephemeris. Colour corrections 
for deriving monochromatic fluxes from the in-band S7 and 
Sll magnitudes can be significant owing to the widths of 
the filters and we used those for a 800-K black body given 
by Lorente et al. (2008, table 4.8.11). There are no near-IR 
data at this photometric phase, so we interpolated J and Ks 
from the ISPI and IRSF observations and show the result- 
ing SED in Fig.m The black body fitted to the difference 
between this and the wind SED has Tbt = 834 ± 25 K with 
r.m.s. residual 0.05 mag. A fit without the Sll data gives 
the same temperature, so we are confident in assigning the 
5*11 fiux to the WR star. Realistic modelling of this SED 
would require a radial dust density distribution to reflect 
the varying rate of dust formation in the preceeding months 
but it is evident from the black-body fits that the character- 
istic temperature of the dust is now significantly cooler than 
that of the phase 0.9 SED, as expected from observations on 
the fading branch of the light curve. 



4 DISCUSSION 

Although the IR light curves need strengthening with fur- 
ther observations, they provide enough evidence for us to 
consider HD 36402 as a variable dust-making CWB like 
WR137. The orbit of WR137 (Lefevre et al. 2005) has pe- 
riastron passage ~ 1.3 ± 0.5 yr before near-IR maximum 
(Williams et al. 2001), which is when the condensation of 
new grains ceased, or at least fell below the rate needed to re- 
plenish the dust being carried away in the stellar wind. This 
puts periastron in WR 137 close to the middle of its 3-yr dust 
formation episode. The dust formation by HD 36402, with 
condensation beginning some time before IR phase 0.7 (cf. 
Fig.[3]l and ending some time near phase 0.0, looks very sim- 
ilar and suggests dust formation formed in a CWB having 
an elliptical orbit with periastron passage near our IR phase 
0.85. The candidate CWB comprises the 08I supergaint in 
an 'outer' orbit with the 'inner', 3.03-d WC4-I-0? binary, 




1 2 5 

Wavelength (ftm) 

Figure 4. SEDs of HD 36402 near photometric phases 0.9 and 
0.2. The dotted line (blue in the on-line version) represents the 
stellar wind spectrum fitted to v and / and determined from that 
of the dust-free system HD 35517. The </> ~ 0.9 SED based on 
the 2001.09 2MASS, 2005.82 IRAC and 2010.41 WISE data is 
marked ffi. A simple model adding a best-fit black body to the 
wind emission is shown with a solid line. The (/) ~ 0.2 SED, based 
on AKARI and interpolated J and Ks photometry is marked •, 
and the corresponding model is drawn with a dashed line (red 
on-line). 



whose orbit was determined by Moffat et al. (1990). They 
found that the absorption lines attributed to the 08I super- 
gaint did not share the motion of the WC4 star, suggesting 
that it was a third, more distant member of the system. We 
re-examined the absorption-line RVs measured by Moffat et 
al., noting that they were observed in three widely spaced 
runs each spanning 3-8 d. We determined the average RV 
from each run to be 283±5, 299±8 and 241±7 km s"^ at IR 
phases 0.17, 0.41 and 0.83 respectively. To test whether these 
represented orbital motion of the 08I star about the 'inner' 
WC4-I-0? binary, we examined the corresponding emission- 
line RVs measured by Moffat et al., adopting their orbital 
elements apart from the 7 velocities as givens, and solving 
for the mean 7 velocity for each run to look for motion of 
the inner binary as a single object. The mean 7 velocities for 
the three runs also showed significant dispersion, but not in 
anti-phase with the absorption-line velocities as one would 
expect from orbital motion, so orbital motion of the 081 star 
has yet to be determined. Further spectroscopic observations 
on a time-scale of years are needed to search for search for 
motion and derive an orbit, but here we will continue on the 
basis that one exists. 

To estimate the properties of the 'outer' (4.7-yr) CWB, 
we use the velocities from the lUE study by Fitzpatrick et 
al. as terminal velocities (woo) and adopt mass- loss rates for 
the WC4 and 08I components of 2.5 x lO"'' and 5.3 x 10"® 
M© y~^ from those of WC4 and 081 stars given by Crowther 
et al. (2002) and Markova et al. (2004) respectively. These 
give a wind-momentum ratio r) = {Mvoo)os/{Mvao)wc — 
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Table 3. X-ray observations of HD 36402 with ROSAT PSPC 
(rp), ROSAT HRI (rh) and XMM-Newton EPIC. The lumi- 
nosities Lx are 'absorbed luminosities', uncorrected for extinc- 
tion; those of the ROSAT observations were derived adopting 
kT = 1.9 keV derived from the XMM-Newton spectral fit. 



Obs ID Date Exp Count rate Lx IR. 

ks ks"-*- erg s"-*- <f> 



source moves through the WC4 wind according to the orbital 
elements - expressions are given by Williams et al. (1990) 
- and is usually greatest near periastron, when the source 
is most deeply buried in the WR stellar wind, or conjunc- 
tion, when the source passes behind the WR star. The non- 
detection of X-ray emission from HD 36402 at phase 0.99 
and relative weakness at phase 0.81 are consistent with this 
picture but modelling the X-ray variations requires knowl- 
ege of the 'outer' orbit. The intrinsic emission measure is 
also expected to vary round the orbit, as 1/D where D is 
the separation of the stars, for an adiabatic WCR (Stevens, 
Blondin & Pollock 1992). This would make the intrinsic 
emission measure somewhat greater at periastron, but with- 
out more spectroscopy or at least measures of the hardness 
ratio around the orbit, we cannot tell at this stage. 

The relatively wide separation of the outer binary im- 
plied by its period should allow non-thermal radio emis- 
sion from the WCR to emerge (cf. Dougherty & Williams 
2000), at least away from periastron, but the distance of the 
LMC would make this a challenging observation with cur- 
rent instrumentation - e.g. the maximum 6-cm flux density 
of WR 140 shifted to the distance of the LMC would be ~ 
30 nJy. 

The demonstration that one of the WR systems in the 
LMC is an episodic dust-maker like those in the Galaxy 
suggests that there are more to be found. Other candidate 
dust-variables include HD 32125 (BAT99-9) and HD 38448 
(BAT99-125). In the case of HD 32125, most of dispersion 
in Ks comes from the DENIS observation, which has a rel- 
atively large uncertainty (a 0.22 mag), so it would be pre- 
mature to consider it a episodic dust-maker and both stars 
need re-examination for variability when photometry from 
the VMC survey becomes available. If episodic dust forma- 
tion is possible in the low metallicity environment of the 
LMC, what about persistent dust formation, like that ob- 
served from some WC8-9 stars observed in the inner re- 
gions of the Galaxy? One candidate amongst the LMC WR 
systems is BAT99-121, which has anomalously red J-Ks 
colours in both the 2MASS and IRSF MC surveys. This 
star is in a difficult field and will be discussed elsewhere. 
In any event, HD 36402 has posed the question: why has 
dust formation not been observed from Galactic WR sys- 
tems outside the solar circle, where the overall metallicity 
is lower? A second question is what is the lowest metallicity 
environment in which WR systems can make and contribute 
dust to the ISM? 
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0.1. Consequently, the distance between the stagnation 
point, where the dynamical pressures of the winds bal- 
ance, and the 08I component will be about one-quarter 
Wv/i'^ + Vv)) that of the separation of the 081 and WC4 
stars. This places the WCR close to the 08 supergaint and 
far from the 'inner' binary, probably accounting for the lack 
of modulation of the dust emission in the WISE data (Sec- 
tion [2TT] above) with the 3.03-d period of the latter. 

Bomans et al. (2002) drew attention to the high X- 
ray luminosity {Lx = 9 x 10^^ erg s"^ of HD 36402 mea- 
sured with XMM-Newton, but noted also the softness of the 
spectral index (P — 2.6), which argued against its being 
caused soley by collision of the stellar winds. X-ray emis- 
sion by HD 36402 was also observed using the ROSAT 
Position Sensitive Proportional Counter (PSPC) and High 
Resolution Imager (HRI) with 0.5-7.0 keV luminosities of 
5.4 ± 1.4 X 10^^ erg s"^ and 3.2 ±0.9 x lO^"* erg s"^ respec- 
tively (Guerrero & Chu 2008). The individual 7? OSAT ob- 
servations were recovered and are given in Table [S] To facili- 
tate comparison, luminosities were derived from the ROSAT 
count rates using web pimms and adopting the absorption 
column Nh = 3.4 x 10^° cm"^ and kT = 1.9 keV derived 
from the XMM-Newton spectral fit, i.e. folding all the vari- 
ability into the 'absorbed luminosities' Lx- 

It is apparent that the X-ray flux observed from HD 
36402 is variable. Most striking are the non-detections in 
the 1992 ROSAT PSCP observations in 1992.19 although 
the integration times were longer than those before and af- 
terwards, in which HD 36402 was detected. The 1992.19 
observations occurred close to zero IR phase, and perhaps 
~ 0.15P after periastron passage in the 4.7-yr binary be- 
lieved to be responsible for the variable dust formation. The 
non-detection of X-ray emission from HD 36402 followed by 
its quick recovery invites comparison with the narrow mini- 
mum in the X-ray light curves of WR 140, which occurs just 
after periastron passage (e.g. Pollock et al. 2005). 

The X-ray fiux observed from a CWB depends on both 
the intrinscic emission measure and the photoionization ex- 
tinction through the stellar winds, especially the heavier 
wind of the WC star. The source is expected to lie close to 
the stagnation point between the WC4-f O? and 08I winds 
and, as seen above, much closer to the 081 star. The cir- 
cumstellar extinction varies as the orbit progresses and the 
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